Temperature-dependent Drude transport in a two-dimensional electron gas 
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We consider transport of dilute two-dimensional electrons, with temperature between Fermi and 
Debye temperatures. In this regime, electrons form a non-degenerate plasma with mobility limited 
by potential disorder. Different kinds of impurities contribute unique signatures to the resulting 
temperature-dependent Drude conductivity, via energy-dependent scattering. This opens up a way 
to characterize sample disorder composition. In particular, neutral impurities cause a slow decrease 
of conductivity with temperature, whereas charged impurities result in conductivity growing as a 
square root of temperature. This observation serves as a precaution for literally interpreting metallic 
or insulating conductivity dependence, as both can be found in a classical metallic system. 

PACS numbers: 72.10.-d 73.40.-c 71.30.4-h 81.05.Uw 
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I. INTRODUCTION 

Electron transport is conventionally understood within 
the Fermi liquid theory framework.^ In a Fermi liquid 
(FL), screening of the offset charge is very efficient due 
to large density of states. As a result, effective disor- 
der potential for quasiparticles is always short-ranged. 
The latter leads to temperature-independent Drude con- 
ductivity, while interaction effects provide corrections to 
the Drude transport in the powers of small parameter 
T/tp, where T is temperature [ks = 1), and ep is the 
Fermi energy. These corrections originate due to scatter- 
ing off the Friedel oscillations,ii^ and due to temperature 
dependence of the random-phase approximation (RPA) 
screening 

What happens when the carrier density n is reduced 
so much that the system becomes non-degenerate, T > 
tpl In this case the Fermi energy is irrelevant, and the 
system is essentially a classical plasma (Fig. [Ij. Such 
a situation may occur in actively studied clean dilute 



heterostructures^'^'i'^'^'i^ where n < 10"'^'^ cm , with 
record densities down to 7 x 10 cm^ -2 (Ref. These 
densities correspond to ep ^ IQ — 100 /iV^ 0.1 — IK. 
Transport in such a system will depend on the strength 
e^/a of electron interactions relative to temperature (here 
a = \l ^fwn is the Wigner-Seitz radius, and the dielectric 
constant k is included into the definition of charge, 
c^/k, for brevity). For strong interactions, e^/a > T, we 
have a strongly-correlated semiclassical electron "liquid" 
whose collective modes are likely to affect transport. ^^'^^ 
Here we consider the two-dimensional (2D) transport 
in the opposite, weakly-interacting regime. 
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We assume that temperature is high enough so that carri- 
ers form a classical weakly-interacting plasma, yet is well 
below the Debye temperature so that the phonon 
contribution to transport can be either neglected or sub- 
tracted in a controlled way. Transport is then dominated 
by the practically unscreened potential disorder. Such a 
situation can become relevant in the cleanest heterostruc- 
tures (e.g. Ref. H), where, for the lowest densities, the 



interaction energy e^/a « 5K is one order of magnitude 
below the Debye temperature. With increasing sample 
quality, the carrier density decreases and the applicability 
range ([1]) widens. Another system where the present ap- 
proach may be applicable is graphene with the substrate- 
induced gap, as described towards the end of the paper. 

We show that in the regime ([T]), the Drude conduc- 
tivity (T(r) becomes strongly temperature-dependent. 
Its temperature dependence originates from the energy- 
dependent impurity scattering cross-section. Remark- 
ably, different kinds of potential impurities (e.g. charged, 
neutral) can now be distinguished by qualitatively dif- 
ferent energy dependence of scattering, yielding unique 
signatures in the resulting <t{T). These signatures could 
be used to characterize notoriously unknown potential 
profile for high-quality 2D samples. 

In particular, for the important example of charged im- 
purities within a 2D layer, we show that the conductivity 
grows as a (X Vt as long as temperature is below a few 
Rydberg of the host material (Ry = TOe^/2?i^ where m 
is the effective carrier mass), crossing over to a oc T for 
T > Ry (Fig. [J). The latter hnear T-dependencei^ is 
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FIG. 1: 2D transport in the presence of charged disorder 
for rs — (e^/a)/eF > 1. FL calculations^'^'^ provide small 
corrections to the Drude transport at T <^ ep. For non- 
degenerate carriers, the perturbative result a oc T^'^ is valid 
for r > Ry , e^/a, while for e^/a < T < Ry we show that 
a oc r^''^. For ef ^5 ^ < e'^/a the system forms a strongly- 
correlated plasma where both Drude^^ and hydrodynamio^^ 
effects can be relevant. For e^/a < ep, i.e. < 1, weakly- 
interacting Fermi gas crosses over to weakly-interacting clas- 
sical plasma at T ~ ep; the result a oc T^/^ then holds for 
EF < r < Ry . 
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thereby practically inobservable for two-dimensional elec- 
tron gases (2DEGs) in GaAs heterostructures, since the 
phonon contribution dominates above 6_d « Ry « 60K. 
Hence, the single-particle explanation of Das Sarma and 
Hwang^"^ of the observed^i^!^'^ conductivity increase with 
temperature does not apply. For the other practical ex- 
ample, the strong neutral impurities, the conductivity is 
shown to decrease with temperature (as described be- 
low). 

As a result, the superficial distinction between a 
"metal" {da/dT < 0) and an "insulator" {da/dT > 0) 
based simply on the sign of the derivative da/dT, does 
not hold - indeed, both behaviors are possible in a clas- 
sical 2D metal ([T]). Of course, a true insulator is char- 
acterized by localized states as T ^ 0, leading to acti- 
vated conductivity dependence. Such a low-T analysis is 
beyond the scope of this work which considers only suffi- 
ciently high temperatures above the onset of localization. 

In what follows, we first obtain the general result ([6|) 
for the T-dependent Drude conductivity in the regime 
([1]) in terms of the energy-dependent impurity transport 
cross-section Atr(e), then we discuss the resulting a{T) 
for different kinds of potential disorder, and, finally, re- 
mark on the systems where one can practically observe 
the temperature-dependent Drude conductivity. 



II. THE DRUDE TRANSPORT 

The kinetic equation in the presence of an external in- 
plane field £ 

vajo = -T-^Sf , e = (2) 
is written in terms of the momentum relaxation rate 



7iiAtr(e)w(e) , Atr = 



^6*^(1 -cos 61). (3) 
du 



Here dK/dO is the differential scattering cross-section and 
Hi is the area density of impurities. The particular en- 
ergy dependence of the scattering rate stems from 
that of the transport cross-section Atr. The isotropic dc 
conductivity follows^"*: 



jde €T^{-defo) 
Jdee{-dJo) 



(4) 



where we assumed energy-independence of the 2D den- 
sity of states in the case of the parabolic band. In the 
classical regime ([T]), quantum interference effects are ir- 
relevant due to strong dephasing. As long as T ^ e^/a, 
one can also neglect electron-electron interactions, such 
that the equilibrium velocity distribution is Maxwellian: 



fo{e{v)) 



,(m-£)/t 



(5) 



Eqs. (gl) and jS]) yield 



a{T) = ao / ^d^ 
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FIG. 2: Temperature-dependent Drude conductivity (|6]) in 
the units of ao = {e'^ /h)n/ni. Solid blue line: charged im- 
purities, Eq. I|lip . together with its asymptotic limits (thin 
dashed lines). Switching from a oc T^''^ to a (x T occurs at 
T ~ 4.5 Ry. Dashed red line: neutral impurities (imperme- 
able disks with radius a = as)- Dash-dotted green line: 50% 
charged and 50% neutral impurities, with the same total n^. 



where energy-dependent wavelength A(e) = 27r/fc(e), 
hk — mv{e). In other words, the temperature depen- 
dence of the Drude conductivity is determined by the 
energy dependence of the transport cross-section in the 
units of wavelength. For simple estimates, Eq. © gives 



cr(T) ~ (To 



A7 



A. 



tr e=T 



(7) 



where At = 2nh/y2mT is the temperature wavelength. 

When multiple kinds of impurities are present, the 
scattering rates add up according to the Matthiessen rule. 
Thus the total transport cross-section entering Eq. ^ 

AUe)=c,A[l\e) + C2A[fie) + ... (8) 

where Cj = rif^ jui is the fraction of impurities of the sort 
j, and Ui — rt'p is the total impurity concentration. 

III. CHARGED IMPURITIES 

For the e^/r potential, the exact 2D differential cross 
section has been found in the seminal 1967 work of Stern 
and Howardi^: 



dbi!^ a tanh ira 
IF ~ 2fcsin2(6l/2) 



a{v) 



e 

fiv 



(9) 



Here 9 is the scattering angle, and the momentum trans- 



fer q 



2hk sin ■ 



The result 



has two distinct 
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limits. For small energies, e <C Ry , the parameter 
7ra » 1, and the cross section is classical (indeed, it is h- 
independent when tanhTra =1). Conversely, for high en- 
ergies (e ^ Ry ), the cross-section ^ with tanhTra ~ Tra 
coincides with the Born approximation. Such a classical- 
to-quantum crossover is specific to 2D, whereas in 3D the 
Rutherford cross-section coincides both with the classical 
result and with the Born approximation.^^ 

The corresponding 2D transport cross-section ([3]) reads 



Ajj, = (2Tra/fc) tanhTra . 



Notably, it is finite, with all scattering angles contribut- 
ing roughly equally. This should be contrasted with 
the well-known logarithically divergent transport cross- 
section for 3D Coulomb plasma-^ ("Landau logarithm"), 
dominated by forward scattering processes. 

The conductivity then follows from Eqs. © and PT!]) : 



i tanh T 



Ry 



(11) 



The conductivity pT|) grows with T (Fig. [5]) since impu- 
rity scattering ([9]) weakens for faster moving carriers. 

The asymptotic behavior of PT|) is a'^ ~ jao^/nTfRy 
at T < Ry , and a"" ~ ^aoT/Ry at T » Ry . Practi- 
cally, the switching between the two limits occurs when 
T fa 4.5 Ry (Fig. ^ the two asymptotes cross). For 
T ^ 5 Ry we agree with Ref. [ij where the Born approx- 
imation TjT^ w rii-K^e^ /fie was utilized in Eq. ^ [corre- 
sponding to the Fermi Golden Rule]. The limit a cx VT 
is novel and relevant in the parameter range ([T]). 



IV. NEUTRAL SCATTERERS 

For any axially symmetric scatterer, the transport 
cross-section is given byi^ 

2 °° 

Atr = sin2((5„,+i - 5^) , (12) 

m— — oo 

where 5m is the scattering phase shift in the channel with 
orbital momentum m. 

Consider an example of strong neutral 2D scatterers 
within the electron layer. Physically, they can originate 
from interface roughness or neutral atomic defects in a 
heterostructure. For sufficiently small ep and T, the car- 
rier's energy may become much smaller than the poten- 
tial barrier which such a potential creates. It is then 
reasonable to model the scattering potential as being in- 
finitely large within a disk of radius a, and zero out- 
side. In this case, the scattering phase shifts tan S„i = 
Jm{ka) / Ymiko) are given in terms of the Bessel functions 
of the first and second kind, leading to 



A" 



8a/3, 



V4+ln^[2/(7fca)] 



/ca > 1 : 
A;a < 1 . 



(13) 



Here ln7 = 0.577... is Euler's constant. Note the anoma- 
lously efficient scattering at wavelengths A = 2Tr/fc ^ a 
exceeding the impurity size: A"j. ~ A/ln^(A/a) ^ a, 
i.e. the scattering cross-section is determined by the car- 
rier wavelength rather than by the impurity size, thereby 
greatly exceeding the "geometric" limit. This is a known 
universal signature of low-energy 2D scatteringJ^ 
The estimate ([7]) yields 



(10) CT"(r) - doxmin 



In^ 



(14) 



The exact conductivity <t^{T) for strong neutral scatters 
calculated numerically using Eqs. (fT2|) and ([6]), is shown 
in Fig. [2] Its asymptotic behavior for small and large 
T agrees with the qualitative estimate ([Tl)) . In order to 
compare with the Coulomb scattering, we took the disk 
radius a — as to be equal to the Bohr radius as — 

/me^, such that Ca ~ Ry as ^ 10 nm for GaAs. 

We also note that weak short-range scatterers yield 
temperature-independent conductivity. Indeed, the dif- 
ferential cross-section dA/dd = |/(^)P in the Born 
approximation /^""(e) = ~mlJ{q)/h'^^/2^ (Ref. [H) 
yields A^""^" cx A for g-independent formfactor l}{q) cor- 
responding to a short-ranged potential U{r). Eq. ^ then 
results ui a — 2TT<7Q{h'^ /mUY— const. 



V. DISORDER SPECTROSCOPY 

In realistic clean low-density samples multiple kinds of 
disorder, e.g. Coulomb impurities and neutral scatter- 
ers, are present. The conductivity ([6]) and ^ can then 
display a fairly complex sample-specific dependence on 
temperature, governed by relative contributions of dif- 
ferent kinds of scatterers. Fig. [2] shows an example with 
c<= = c" = 0.5. 

Qualitatively different T-dependences (fTTj) and ^4]) 
present a natural way to characterize disorder in clean 
2D samples. For that one needs to operate at very low 
carrier densities n ^ 10^ cm~^, when a temperature win- 
dow ([1]) opens up. Fitting the conductivity (with the 
phonon contribution subtracted) to the result ^ and 

([5]) will yield the disorder composition {np^}. This 
way, the T-dependent transport can serve as the dis- 
order spectroscopy. The connection with spectroscopy 
is not accidental: Formally, the conductivity is propor- 
tional to Laplace transform J^d€e~^'^ip{e) of the quan- 
tity ^(e) = Vi/Atr(e), where /? = 1/T. 



VI. CHARGED DISORDER IN GaAs 
HETEROSTRUCTURES 

The result pT|) based on the exact cross-section 
predicts a novel a oc Vt conductivity dependence, 
characteristic of the classical limit of scattering 
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The latter can be relevant for transport in clean di- 
lute heterostructures4i^i^i2ii2ii2iii^ So far, the observed con- 
ductivity grows approximately linearly with temperature 
around T ~ 1 From the present analysis, the 

single-particle explanation-'^'^ for this observation based 
on the Born scattering cannot hold for GaAs, since, ac- 
cording to Fig. [21 the crossover to the Born regime would 
occur at T « 300 K which is practically inaccessible due 
to the dominant phonon scattering j^^i^^ The apparent 
discrepancy between the present single-particle theory 
yielding a cx VT, and the experiments^'^'i strongly in- 
dicates the predominance of collective effects in trans- 
port. This is not surprising, since typical Coulomb en- 
ergy e^/a ~ 20K [corresponding to n = 1 x lO^'^cm^^], 
while the measurements were done for at least order-of- 
magnitude lower temperatures, in which case using the 
Maxwell distribution (O in Eq. (|4]) is unjustified from 
the outset. For the lower densities, n ^ lO^cm"^, the 
present approach may apply, as long as the phonon con- 
tribution is controllably subtracted in the range ([T|). 

Can the linear (RPA) screening affect the temperature 
dependence (jlip . and in particular, the crossover tem- 
perature T w 5 Ry ? Below we argue that screening will 
only weaken the dependence cr{T), and cannot lead to 
cr(r) cx T at low temperature. 

Physically, screening changes the shape of the impu- 
rity potential in the following way: It fully preserves 
the strength of the e^/r potential for distances r < 
shorter than the screening length, and cuts off the 1/r 
behavior for r > a,, where a, = „ -^o = ^-rr- fin the 
Fourier space, iire^ /k 27re^/(fc + a^^)]. The linear 
(RPA) screening is a mean-field effect, valid when the 
density fluctuations within the screening volume are 
small, fulfilled under the condition naj. ^ 1 equivalent 
to T ^ e^/a, compatible with the limit ((!]). This has the 
following consequences: (i) for < T < e^/a relying 
on the RPA screening is unjustified. The single-particle 
transport calculation based on the Maxwell distribution 
([5]) is also unjustified. Thus the approach^'^ of Das Sarma 
and Hwang does not apply to the experiments^'^'i even 
if the authors were to use the correct scattering cross- 
section, (ii) For T ^ e^/a, screening becomes asymp- 
totically irrelevant for the Drude transport. Indeed, con- 
sider the region r < where the electron "feels" the un- 
screened impurity potential. Upon entering this region, 
its typical kinetic energy greatly exceeds the Coulomb 
field, T ^ e'^/as- Thus the scattering phase shifts yield- 
ing the cross-section © have parametrically large room 
to accumulate between e^/T <C r <C a^, leading to its 
nonperturbative limit. Moreover, the residual screening 
(truncation of the potential for r > a^) would further 
weaken the (j{T) dependence, since, according to the 
above calculation [cf. Eqs. (O and (O], the conductivity 



due to short-range disorder decreases with temperature. 
Thus the initial a cx Vt dependence would only weaken 
when the residual screening is taken into account. 

As a result, the explanationi^ suggested for the 
apparent linear growth of the conductivity with 
temperature^i^i^ii does not apply; the observed linear 
(and, generally, power-law*) T-dependence of the low- 
temperature conductivity remains an exciting unresolved 
problem. 



VII. GRAPHENE WITH CHARGED DISORDER 

The nonrelativistic scattering considered above can be 
applied to graphene samples where Dirac mass m = 
^/vp can originate e.g. from symmetry- breaking be- 
tween sublattices, such that gap A ~ 10 — 100 meV j^°i^^ 
Half-filled band corresponds to chemical potential fj, = 
—A counted from the bottom of the "parabolic" band. 
The graphene electron system is nonrelativistic and non- 
degenerate as long as r ^ A, since ep = Te~^^'^ <C T 
[Eq. ([5|)]. When electron interactions (controlled by di- 
electric environment) are weak, ag = e'^/hvp ^ 1 where 
VF 10^ m/s, the effective Rydberg Ry^ = a'^gA/2 < A. 
Hence, cf. Fig[l] the conductivity a cx 

for T < Ry 

and (7 cx T for Ry^ ^ T ^ A. For strong interactions, 
ag ^ 1, Ry^ ^ A and the regime a ocT never plays out. 
For T ^ A the system becomes relativistic, the cross- 
section scales as the wavelength^ and the T-dependence 
of the Drude conductivity cr cx comes solely from that 
of carrier density n (xT^ , Ref . [2^ 



VIII. SUMMARY 

To conclude, we considered temperature-dependent 
Drude transport in non-degenerate 2D electron systems. 
The Drude conductivity due to charged disorder behaves 
classically, a cx T^/^ for temperatures below a few Ry- 
dberg, while neutral disorder results in decreasing cr(T). 
These signatures can be utilized in determining disorder 
content of clean 2D samples in the limit (H]) . The decrease 
of the conductivity while reducing temperature does not 
necessarily signify a transition to an insulating state. 
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